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Abstract. We show how, starting from QCD Lagrangian, a low-energy limit can be derived that is 
(NT given by a non-local Nambu-Jona-Lasinio model with the form factor in close agreement with the 

Q_- one expected from an instanton liquid model but with all the parameters properly fixed. In the mean 

QJ ' field approximation, the agreement with experimental data is excellent. 
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Knowing classical solutions to field theories can help to approach the corresponding 
quantum problem. For nonlinear theories these are often difficult to find and one is not 
able to manage a theory other than using some non-perturbative techniques as lattice 
computations or, in a lucky situation, with a small perturbation theory. In this paper we 
aim to show as a class of nonlinear field theories is amenable to analytical treatment, 
providing a class of classical solutions are known to be used as a starting point for 
perturbation theory in the infrared limit of quantum field theory. One of them is Yang- 
Mills theory and so, one has a technique to study this theory in the low-energy limit. 

A low-energy limit for Yang-Mills theory implies that one is able to manage quantum 
chromodynamics (QCD) in the same limit. Indeed, if one knows the propagator of the 
gluonic field in this limit, the behavior of QCD at lower momenta can be obtained [1,2]. 
(Sj | The model that one recovers in such a case is a non-local Nambu-Jona-Lasinio model, 

as already discussed in [3] but with the form factor precisely fixed by QCD. This form 
factor agrees very well with the one obtained through an instanton liquid [4] showing 
that, in the ground state of the theory, instantons could play a fundamental role. It is 
important to note that there exists another approach showing this same link between 
QCD and Nambu-Jona-Lasinio model [5]. 

In the following we show how, starting from a set of classical solutions, the low- 
energy limit of QCD can be recovered yielding a treatable Nambu-Jona-Lasinio model. 
A prototypical field theory is a massless scalar field with a simple non-linearity given 

by 

□0 + A0 3 = ./. (1) 



The corresponding homogeneous equation admits an exact solution § = ji (2/ A) 4 sn(p • 
x+6,i) being sn an elliptic Jacobi function and /! and 6 two integration constants. In 
order for this solutions to hold, the following dispersion relation applies p 2 = /i 2 a/A/2. 
We recognize here a free massive solution notwithstanding we started from a massless 
theory. We see that mass arises from the nonlinearities when A is taken to be finite rather 
than going to zero and so, standard perturbation theory just fails to recover it. Our aim is 



to solve the equation (1) in the limit X — > °° and then we consider an approach devised 
in the '80s [6] taking § as a functional of j and expanding in powers of it. The argument 
is consistent provided one puts 

0(x) = ju J d 4 x'A(x-x')j(x') + fl 2 X J d 4 x'd 4 x"A(x-x')[A(x'-x")] 3 j(x!) + 0(j(x) 3 ) 

(2) 

being A(x — xf) a solution to the nonlinear equation OA(x — x 1 ) + X [A(x — xf)] = 
H~ 1 8 4 (x — x 1 ). So, this current expansion is meaningful as a strong coupling expan- 
sion. The interesting result here is the linear term of the Green function in the current 
expansion. When applied to a quantum field theory this will provide a Gaussian gen- 
erating functional marking a trivial theory. So, we recognize that this set of classical 
solutions yields, at the leading order, a trivial theory. Such a theory is mathematically 
manageable notwithstanding the strong nonlinearity of the theory we started from. This 
program can be accomplished if we know how to get the Green function and this is eas- 
ily seen starting from the d = 1 + theory. In this case we have to solve the equation 
d 2 A (t -t') + X [A (t -t')} 3 = /l 2 8(t - 1') and this is easily obtained [7]. Then, after a 
Lorentz boost or a gradient expansion, the full propagator is given by 

^ % 2 (-l) n e -( n+ b a 1 

Mp) = z>+ D m \ ^ eHMl} , p i_ mi+ie O) 

being m n = (2n+ \){%/2K{i)) (X/2)* /i and K(i) « 1.3111028777 an elliptic integral, 
consistently with the idea of a strong coupling expansion. Then, we would like to apply 
all this machinery to Yang-Mills theory but in order to show this we need to have a set 
of classical solutions to work with also in this case. This set of solutions would grant 
a trivial infrared fixed point also for this theory. Such solutions exist. This can be seen 
starting from the equations of motion 

and assuming again a current expansion. We note that the homogeneous equations 
can be solved by setting A c ^(x) = r\^{x) being r\ c ^ a set of constants. In this case 
we need to fix the gauge and we assume the Lorenz gauge being this equivalent to 
the Landau gauge in quantum field theory. So, the homogeneous equations collapse to 
d^d^ty + Ng 2 ^> 3 = —j<f, and we have turned back to the previous scalar field theory (this 
is no more true for other gauges where the correspondence is just an asymptotic one 
[8]). In this way, the gluon propagator in the Landau gauge is straightforwardly obtained 
from eq.(3) setting X = Ng 2 and with a factor 8 a h {f)iiv — Py.Pv/ P 2 ) ■ 

We started with a fundamental hypothesis to support analytical developments, that is 
that Yang-Mills theory could admit a trivial infrared fixed point for the running coupling. 
Indeed, lattice computations strongly support this view as shown by the German group 
[9] from lattice at 64 4 and 80 4 with jS = 5.7 where the running coupling is seen to 
go to zero as momenta lower. A similar result was obtained by the French group 
with a different definition of the infrared running coupling [10]. They show a perfect 
consistency with an instanton liquid model in agreement with the scenario we are 
depicting here. 



Moving to quantum field theory, the generating functional for the scalar field can 
be managed by rescaling the space-time coordinates as x — > Vhx and with a strong 
coupling expansion = E,7=o^ Then, at the leading order we will have to solve 
the equation D0o + A0q = j that we now know how to manage. Then, the leading order 
is just a Gaussian generating functional with the propagator given by eq.(3) when use 
is made of the approximation in eq.(2), next-to-leading order can be also computed. We 
arrive at the fundamental result that the massless scalar field theory in four dimensions 

is infrared trivial [11]. Mass spectrum is given by m n = (2n + \ ){n/2K{i)) (A/2) 5 ji 
as expected, representing free particles with a superimposed spectrum of a harmonic 
oscillator. Now, turning the attention to Yang-Mills generating functional we realize that 
it takes the simple Gaussian form 



Zo[j] = /Vexp 



[ - J MMj^^nf^ -*W) 



(5) 



once we use the current expansion A a ^ = A J d A x'D a ^ v (x — x') j bv (x') +0(1/ VNg) + 

0(f) and the propagator Df v (p) = 8 ab (tj mv - ^\ A(p) being A(p) given by eq.(3). 

The spectrum in this case is that of free massive glueballs with a superimposed spectrum 
of a harmonic oscillator. This is consistent with our initial observation of a trivial infrared 
fixed point for the running coupling. Similarly, considering the ghost field, applying 
our approximation on the gauge field through instanton solutions, this decouples at the 
leading order producing a free particle propagator for a massless field. These properties 
of the quantum Yang-Mills field describe the so-called "decoupling solution" [12, 13, 
14]. This solution is the one recovered in lattice computations [15, 16, 17] and agrees 
perfectly well with the numerical solution of Dyson-Schwinger equation as given in [12] 
(see fig. 1). 
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FIGURE 1. Comparison of our propagators with numerical solution of Dyson-Schwinger equations as 
in [12]. Running coupling is also shown. 

One can apply this low-energy behavior of Yang-Mills theory to the full QCD gen- 
erating functional, with the propagator the one we just obtained, yielding for the action 
entering into it 



S = 



d 4 x 



(6) 



where the a field arises from the gluon propagator in the Gaussian generating functional 
of the Yang-Mills action [1] and is the contribution from the mass gap of the theory, 
being mo = (n/2K(i))\f& and c is the string tension (440 MeV) 2 ). For the quark 
fields one has 



(7) 



S i = H d 4 xq(x) 

q J 

- g 2 J A'A(x-xOi;^(x)y A'(xOyrM?'(^?W + o +0(j 3 ). 

Now, we are able to recover the non-local Nambu-Jona-Lasinio model given in [3] but 
directly from QCD provided the form factor is 



Sf(p) 



1 



g 2 Mp) 



l 



B,, 



2 6 ^ p 2 -(2n + l) 2 (n/2K(i)) 2 d + ie 2 



V(p) (8) 



being B n obtained from eq.(3), ^(0) = 1 and 2^(0) = G the standard Nambu-Jona- 
Lasinio coupling, fixing in this way the value of G through the gluon propagator. We 
compare this form factor both with the one guessed in [3] and the one from an instanton 
liquid [4] in fig. 2. The result is strikingly good for the latter showing how consistently 
our technique represents Yang-Mills theory through instantons. 




FIGURE 2. Comparison of our form factor with those provided in [3] and [4] as shown in [1]. 



Finally, we take the local Nambu-Jona-Lasinio limit in the non-local model and 
compute some observables through equations given in [18]. The agreement is absolutely 
excellent as seen in tab. 1, notwithstanding the rough approximation, for A = 1 183 MeV 
(4d cut-off), 6 = (440 MeV) 2 , m q = 4.1 MeV and g = 1.52(a, « 0.18). Then, G « 
0.7854(g 2 /^) « 4.9 GeV 2 that gives G eff w 2.6 GeV 2 and so G eff A 2 w 3.6 being 
< G the Nambu-Jona-Lasinio model modified by the presence of the 



J eff 



mass gap arising after bosonization of the action (7) as in [3]. We use mean field formulas 
as given in Klevansky [18] and refer to PDG for experimental values. 

We have shown how, starting from QCD Lagrangian, it is possible to derive a non- 
local Nambu-Jona-Lasinio model with all parameters properly fixed. We were able to 
compute some meaningful observables in close agreement with experimental data even 
if in a mean field approximation. 



TABLE 1. Comparison with experimental values in the contact approxi- 
mation [18] 



Observ. 


Exp. 


Theor. 


Error 




139.57018±0.00035 MeV 


139.7 MeV 


0.09% 


f* 


130.41±0.03±0.20MeV 


131.5 MeV 


0.9% 


-(mm) 3 = -(dd)? 


230±10MeV (sum rules) 


274 MeV 


16% 


gnqq 
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